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Abstract 

The volume dependence of the octet baryon masses and relations among them are explored with 
Lattice QCD. Calculations are performed with nj = 2 + 1 clover fermion discretization in four 
lattice volumes, with spatial extent L ~ 2.0, 2.5, 3.0 and 3.9 fm, with an anisotropic lattice spacing 
of b s ~ 0.123 fm in the spatial direction, and bt = b s /3.5 in the time direction, and at a pion mass of 
m-Tr ~ 390 MeV. The typical precision of the ground-state baryon mass determination is < 0.2%, 
enabling a precise exploration of the volume dependence of the masses, the Gell-Mann-Okubo 
mass relation, and of other mass combinations. A comparison of the volume dependence with the 
predictions of heavy baryon chiral perturbation theory is performed in both the SU(2) L (g) SU(2) R 
and SU(3) £ (g) SU(3) R expansions. Predictions of the three-flavor expansion for the hadron masses 
are found to describe the observed volume dependences reasonably well. Further, the AA^7r axial 
coupling constant is extracted from the volume dependence of the nucleon mass in the two-flavor 
expansion, with only small modifications in the three-flavor expansion from the inclusion of kaons 
and rj's. At a given value of m n L, the finite- volume contributions to the nucleon mass are predicted 
to be significantly smaller at ~ 140 MeV than at m n ~ 390 MeV due to a coefficient that scales 
as ~ m 3 . This is relevant for the design of future ensembles of lattice gauge- field configurations. 
Finally, the volume dependence of the pion and kaon masses are analyzed with two-flavor and 
three-flavor chiral perturbation theory. 
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I. INTRODUCTION 



The calculation of the properties and interactions of light nuclei is a major goal of Lattice 
QCD. While Lattice QCD calculations at the physical light-quark masses, including strong 
isospin breaking and electroweak interactions, are a number of years in the future, precision 
calculations of hadron masses are being performed today in the isospin limit and without 
electroweak interactions over a range of light-quark masses. The masses of the baryons and 
nuclei are in the GeV energy regime, but the typical excitation energies and binding energies 
found in light nuclei are in the MeV energy regime. This hierarchy presents a significant 
challenge for Lattice QCD calculations as correlation functions must be determined with 
exceptionally high precision in order to obtain statistically significant energy differences 
that yield nuclear excitation and binding energies. 

A source of systematic uncertainty in the extraction of scattering parameters and nuclear 
binding energies is the volume dependence of the hadron masses themselves. Given that 
the deuteron binding energy is Bp ~ 2.2 MeV, an accurate Lattice QCD calculation of 
this energy will require that the nucleon mass be known to a precision of AMn <C 1 MeV. 
This includes the contribution from the finite lattice volume. Further, the exponential 
volume corrections to Liischer's eigenvalue relation [THS] are also required to be small [HE]. 
In our recent calculation of the H-dibaryon binding energy [6], the volume dependence of 
the A-baryon mass was presented, and it was concluded that the standard rule-of-thumb, 
m-n-L^ 2tt, is in fact necessary at a pion mass of m n ~ 390 MeV in order for the A finite- 
volume mass shift to be much smaller than the observed binding energy. 

In this work, which is a continuation of our high statistics Lattice QCD explorations [3-0] , 
we present results for the volume dependence of the masses of the baryons in the lowest- 
lying SU(3)-flavor octet, and of relations among them, calculated with four ensembles of 
n/ = 2 + 1 anisotropic clover gauge- field configurations at a pion mass of m n ~ 390 MeV 
with a spatial lattice spacing of b s ~ 0.123 fm, an anisotropy of £ = 3.5 and with cubic 
volumes of spatial extent L ~ 2.0, 2.5, 3.0 and 3.9 fm. The volume dependence of the pion 
and kaon masses are also determined. Having the multiple lattice volumes with all of the 
other parameters fixed is critical to fully understanding the volume dependence of the hadron 
masses and other quantities. In particular, lattice-spacing artifacts, which enter at 0{b 2 s ), 
are the same in all four ensembles to very high precision. The results of the Lattice QCD 
calculations are compared with the expectations from next-to-leading-order (NLO) chiral 
perturbation theory (%PT) and heavy-baryon chiral perturbation theory (HB^PT) with 
two and three flavors of active quarks l . While it is interesting to compare the calculated 
volume dependences with the corresponding expectations from low-energy effective field 
theories (EFT), perhaps the most important reason for such a study is in order to plan 
for the future production of ensembles of lattice gauge-field configurations. An interesting 
result of this study is that while the octet baryons experience significant and quantifiable 
finite- volume corrections, the pion and the kaon, whose masses are determined at the 0.1% 
level, show effectively no finite-volume effects. That is, the pions and kaons are in the 
infinite-volume regime for the pion mass and the range of volumes that we explore. These 
results are consistent with the expectations derived from %PT and HB^PT. 



The leading volume dependence of the hadron masses arises at NLO in the chiral expansion. The leading 
order (LO) contributions to hadron masses from local operators make vanishing contributions to the 
volume dependence. 
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The first realistic attempt to determine the coefficients of counterterms in the chiral 
Lagrangian from the volume dependence of the nucleon mass was performed in Ref. [TO] . 
The coefficients in the SU(2) L <g> SU(2) R chiral Lagrangian without dynamical A's were con- 
strained by the results of rif = 2 Lattice QCD calculations using the clover discretization 
with m n > 550 MeV and with lattices of spatial extent L< 2.2 fm. The precision of these 
calculations was substantially lower than in the present work, nonetheless, nontrivial con- 
straints were found on the values of coefficients in the chiral Lagrangian at NNLO in the 
expansion by using phenomenologically determined values for the NLO constants. These 
constraints should be viewed only as a demonstration of the method, given the large pion 
masses. 

This paper is organized as follows. In section [TTJ we formulate finite-volume correction 
formulas for octet-baryon masses to NLO in SU(2) L ®SU(2) R and SU(3) L ®SU(3) fi HBxPT. 
Section III gives a concise description of the specific Lattice QCD calculations that are 
used in the present finite- volume study. In section |IV[ we analyze the octet-baryon finite- 
volume effects, first (in subsection IV A) using a simple, intuitive description, and then (in 
subsection IV B) using HByPT. In subsection IV C various combinations of baryon masses 
are likewise analyzed. In section [V] we consider the finite-volume dependence of the pion 
and kaon masses and in section [VI] we conclude. 



II. FINITE- VOLUME CHIRAL PERTURBATION THEORY 

A. The Nucleon in SU(2) L ® SU(2) iJ HB^PT 

As the Goldstone bosons are the lightest hadrons, is the appropriate tool to develop 
systematic expansions to describe finite- volume effects [TTHT3] . The crucial observation is 
that if the hadronic system is in a sufficiently large volume, then the infinite- volume chiral 
Lagrangian can be used to calculate finite-volume corrections, with no further operators 
required 2 . At NLO in HB%PT, the finite-volume corrections to the nucleon mass arise from 
one-loop self-energy diagrams with nucleon and A intermediate states and, in the limit of 



The finite- volume corrections are related to forward scattering amplitudes 14-16 . For instance |15j . 



M N (m w ,L) -M N (m w , oo) = M N V ri ^— r \ 2*e*& N e~^ m - L 

n^O 1 1 W 



+ (v) P -\n\m^L \/ 1 +y 2 



dy D + (y) e 



+ 0{e~ ML ) , 



(1) 



where — m 7r /(2Afjv), M > y / 372m 7r and D + (y) = M/v D + (im 7r y, 0), which is related to forward 7T./V 
scattering via 



1 



T(n a (q) + N(p)^7r a (q')+N{p')) = S aa ,T+ + -[r a ,,r a } T 



D ± (v, t) 



AM 



N 



B ± {u 1 t) 



(2) 



The strong coupling between the nucleon and pion is g^N, which is related to the axial coupling constant 
via the chiral expansion g^N — 9aM^-\/2/ + ... where f„ ~ 132 MeV. Evaluating eq. ([IJ at NLO in 
HB^PT recovers the perturbative result of eq. (pi). 
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exact isospin symmetry, are given by [TOl [T7] 3 

SMn = MN(m n , L) — M]y(m n , oo) ; 

9mlg 2 A „(fv), n m l9 2 A N^ Ufv) 



F$ v > (m^L) + ^ v , AanL) j (3) 



where the loop functions are defined to be 

n^O 



= rdw/3(w,V) /3(w^)K (\n\x/3(w,^ 

Sir In x x \ x 

n^O Jl) 



Ki{\n 






n 


a; 



(5) 



where f3(w,z) = \/w 2 + 2zw + 1, Aatv = Ma — M^v, and the K n (z) are modified Bessel 
functions. In the limit A — > 0, F^' V ^(m 1T L,AL) — > F^f v \m 7T L). For asymptotically large 
lattice volumes, it is expected that only the leading contributions in the sums in eq. Q and 
eq. ([5]) will be relevant. 

B. The Hyperons in SU(2) L ® SXJ(2) R HB^PT 

It is straightforward to compute the finite- volume corrections to the hyperon masses at NLO 
in SU(2) L eg) SU(2) fl HB%PT. One finds, in a generalization of eq. ([3]), that 

SM A = 3 -4^Fi FV \ m7T L,Az A L) + 3 %^ Ff^(m 7r L,A^ A L) ; (6) 

+ %#^r ) (-^A s , E L); (7) 

where Aab = Mb — Ma- For definitions of the various axial couplings in terms of chiral 
Lagrangian operators, see Ref. [I8j. In this formulation of hyperon finite- volume effects, the 
contributions from kaon and rj loops are in the coefficients of local operators, and therefore 
do not contribute until higher orders in the chiral expansion (as the finite- volume effects arise 
from pion loops). This implies that for quark masses such that m v < (but not m n <C 
tuk) there is a region where the two-flavor chiral expansion of infinite- volume quantities will 
converge but finite-volume effects will not be accounted for systematically in the two-flavor 
expansion. 



3 The substitutions w — > A/m^, z — > A/m T , P(w, |) — > /3aA?V, V ~> AL, x — > m n L, and F^ V \x 1 y) 
— /C(A)/(37rm 3 r ) recover the expressions given in Ref. [17] , 
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C. The Baryon Octet in SU(3) L <g> S\J(3) R HBxPT 

In addition to the contributions to the volume dependence from higher orders in the 
SU(2) L sg) 811(2)^ chiral expansion, there are contributions from quantum fluctuations of 
the nucleon into strange hadrons, and of the hyperons into other members of the baryon 
octet. For instance, in addition to the nN and 7rA intermediate states that give finite 
volume contributions to the nucleon mass, intermediate states such as AK or Nrj also con- 
tribute. Such fluctuations give rise to a volume dependence that scales as ~ m 2 K e~ mK ' L /L 
or ~ m 2 e~ mr > L /L. In the Lattice calculations presented in this paper, m^/m^ ~ 1.4, and 
as a result, such contributions are naively expected to be of the same magnitude as the 
~ m 2 e~^ m7rL I L contributions, which are suppressed compared with the leading volume 
dependence. As the axial-couplings between the nucleons and pions are of order one, and 
the couplings to strange intermediate states are generically small, such strange contributions 
are expected to be small. The explicit contributions to the octet-baryon mass shifts, written 
in terms of the SU(3)-symmetric axial couplings D, F and C [l9j are for the nucleon, 



5M™> = (D-FY F£ V) (m K L,Ax N L) + (D + F% V) (m K L, A AN L) 

+ (D-3Ff £j- 2 Ff y) KL,0) + C 2 ^| F A FV \m K L,A^ N L) , (9) 
and for the hyperons, 

5M[ K ^ = {D + ZFfp^ F A FV \m K L,A NA L) + (D - 3F) 2 -^ Ff V \m K L, A HA L) 
+ D*^j_ F (^) K L, ) + C 2 ^ F A FV \m K L,A B * A L) , (10) 



SM^ = {D-Ff^Fi FV \ mK L,A N ^L) + (D + Ff^ F A FV \m K L, A HS L) 
+ D 2 ^Fi FV \m v L,0) + C^Fi FV \m v L,A^L) 
+ C 2 !^| Fl FV \m K L,A A xL) + C 2 -^- F A FV \m K L, A S ^L) , (11) 

8M^ = (D + F)^Fi FV \ mK L,A^L) + (D - 3F) 2 -% F A FV \m K L, A AS L) 
+ P + SFf-^Ff^K^O) + C^ 2 Fi FV \m v L,A^L) 

+ C 2 v% F A FV \m K L,A^L) + C 2 ^%- F A F V \m K L, A UE L) . (12) 

In the limit of exact SU(3) symmetry, the SU(2) axial couplings introduced above are related 
to the SU(3) couplings via: 



g A = D + F , g A Nn = C , g A ^ = 2D , g> E * A7r = C/y/2 , 

<fes = 2F , g^ n = C/V3 , 9ss = D-F , g s * Sn = C/V3 . (13) 
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Adding the finite-volume modifications in eqs. (J9|-(12) to those in eq. ^ and eqs. (Jsj) - (J8j) 
gives the full NLO SU(3) L <g> SU(3) R HB%PT finite-volume contributions to the baryons in 
the lowest-lying octet. 



III. DETAILS OF THE LATTICE QCD CALCULATION 

Anisotropic gauge-field configurations have proven useful for the study of hadronic spec- 
troscopy [20H23] . and, as the calculations required for studying multi-hadron systems rely 
heavily on spectroscopy, we have put considerable effort into calculations using ensembles 
of gauge fields with clover-improved Wilson fermion actions with anisotropic lattice spacing 
that have been generated by the Hadron Spectrum Collaboration (HSC). In particular, the 
rif = 2 + 1 flavor anisotropic clover Wilson action [24"l 125] with stout-link smearing [26] of 
the spatial gauge fields in the fermion action with a smearing weight of p = 0.14 and n p = 2 
has been used. The gauge fields entering the fermion action are not smeared in the time 
direction, thus preserving the ultra-locality of the action in the time direction. Further, a 
tree-level tadpole-improved Symanzik gauge action without a 1 x 2 rectangle in the time 
direction is used. 

The present calculations are performed on four ensembles of gauge-field configurations 
with L 3 x T of 16 3 x 128, 20 3 x 128, 24 3 x 128 and 32 3 x 256 lattice sites, with a renormalized 
anisotropy £ = b s /b t = 3.5 where b s and b t are the spatial and temporal lattice spacings, 
respectively. The spatial lattice spacing of each ensemble is b s = 0.1227 ± 0.0008 fm [22] 
giving spatial lattice extents of L ~ 2.0, 2.5, 3.0 and 3.9 fm respectively. The same input 
light-quark mass parameters, b t mi = —0.0840 and b t m s = —0.0743, are used in the produc- 
tion of each ensemble, giving a pion mass of m n ~ 390 MeV. The relevant quantities to 
assign to each ensemble that determine the impact of the finite lattice volume are m^L and 
m n T, which for the four ensembles are m n L ~ 3.86,4.82,5.79 and 7.71 respectively, and 
m^T ~ 8.82, 8.82, 8.82 and 17.64 respectively. 

Multiple light-quark propagators were calculated on each configuration in the four en- 
sembles. The source locations were chosen randomly in an effort to minimize correlations 
among propagators. On the {16 3 x 128, 20 3 x 128, 24 3 x 128, 32 3 x 256} ensembles, an aver- 
age of {224, 364, 180, 174} propagators were calculated on each of {2001, 1195, 2215, 774} 
gauge-field configurations, to give a total number of ~ {4.5, 4.3, 3.9, 1.3} x 10 5 light-quark 
propagators, respectively. The EMP's for the nucleon on all of the time slices of each lattice 
ensemble are shown in fig. [TJ They provide an indication of the precision of the present 
calculations. 

IV. THE VOLUME DEPENDENCE OF THE BARYON MASSES 

The baryon masses that were extracted from the Lattice QCD calculations in the four 
different lattice volumes are given in table |T] A detailed discussion of the fitting methods 
used in the analysis of the correlation functions is given in Refs. [7HS1 EZ] ■ 

A. A Simplistic Analysis 

It is useful to begin the analysis of the volume dependence of the baryon masses by per- 
forming a straightforward, but only partly motivated, fit to results of the Lattice QCD 
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t / b ( t / b, 

FIG. 1: The nucleoli EMP's obtained in the four lattice volumes. Note that the temporal extent 
of the 32 3 x 256 ensemble is twice that of the other three ensembles. 



TABLE I: Results from the Lattice QCD calculations in the four lattice volumes. 



L 3 x T 


16 3 x 128 


20 3 x 128 


24 3 x 128 


32 3 x 256 


L (fm) 
m n L 

1 -m-xL 

1 -ITIttL 


~ 2.0 
3.888(20)(01) 

~0.0205 
~5.3 x 1(T 3 
~2.7 x 1(T 3 
8.89(16)(01) 
~1.38 x 10" 4 


~2.5 
4.8552(84)(35) 
~0.0078 
~1.6 x 10~ 3 
~7.4 x 10~ 4 
8.878(54)(22) 
~1.39 x 1(T 4 


~3.0 
5.799(16)(04) 

~0.0030 
~5.2 x 10~ 4 
~2.2 x 10" 4 
8.836(85)(02) 
~1.45 x 10" 4 


~3.9 
7.7347(74)(91) 
~0.00044 
~5.7 x 10" 5 
~2.1 x 10" 5 
17.679(59)(73) 
~2.10 x 10~ 8 


M N (t.l.u.) 
M A (t.l.u.) 
M s (t.l.u.) 
M s (t.l.u.) 


0.21004(44)(85) 
0.22446(45)(78) 
0.22861(38)(67) 
0.24192(38)(63) 


0.20682(34) (45) 
0.22246(27)(38) 
0.22752(32)(43) 
0.24101(27)(38) 


0.20463(27)(36) 
0.22074(20)(42) 
0.22791 (24) (31) 
0.23975(20)(32) 


0.20457(25) (38) 
0.22054(23)(31) 
0.22726(24)(43) 
0.23974(17)(31) 



calculations given in table [Tj As shown previously, the volume dependence of the mass of a 
given baryon can be calculated order- by-order in HB%PT. The formally-leading contribution 
to the volume dependence of the mass of an octet baryon results from a one-loop diagram 
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involving a pion and an octet baryon (ignoring for the moment the contribution from decu- 
plet baryons). These one-loop contributions give rise to a volume dependence of the form 
given in eq. (4), F^ v \m n L). In obtaining this result, it is assumed that m^L is large, but 
significantly smaller than m^I where mx is the mass of other mesons, such as the kaon and 
the r], i.e. tuk , Tn^ 3> m n . In the very-large volume limit, the finite-volume contributions 
are dominated by the first term in eq. (|4]). As such, it is useful, as a preliminary analysis, 
to fit a function of the form 



.\/; ; x ) r 



(V) e 

B 



(14) 



to the results of the Lattice QCD calculations given in table [IJ One should view the parameter 
as providing an estimate of the strength of the axial coupling between the pion and the 
baryon. It should be stressed that the higher-order terms, beginning with terms of order 
mnL /(m^L), give a non-negligible contribution (relative to the uncertainties) in the 



16 x 128 lattice volume, and a fit using the full function in eq. (|4 



values of compared to those determined in the fits to eq. (|14|). The fits to each of the 



leads to slightly reduced 



baryon masses are shown, along with the results of the Lattice QCD calculations, in fig. [2j 
The same vertical scale (but different intervals) has been used in the plots in fig. [2] in order 
for the reader to easily determine the relative size of the volume dependence of each of the 
masses. The values of the infinite- volume masses, Mg°\ and the coefficients of the leading 
volume dependences, Cg, are presented in table II The nucleon is found to have the largest 



TABLE II: The results of linear fits, of the form given in eq. (14), to the volume dependence of 



the baryon masses. M^°^ is the infinite-volume extrapolation of the baryon mass and is the 
coefficient of e~ mnL / (m^L). The first uncertainty is statistical, the second is the fitting systematic, 
and the third (where appropriate) is due to scale setting. 



Hadron 


(t.l.u.) 


M B °° ] (MeV) 


4 y) (t.i.u.) 


cg° (MeV) 


M N 
M A 
M E 
M s 


0.20427(17)(19) 
0.22053(15)(21) 
0.22744(17)(22) 
0.23972(13)(18) 


1149.8(1.0)(1.1)(7.5) 
1241.2(0.9)(1.1)(8.1) 
1280.3(1.0)(1.1)(8.3) 
1349.4(0.8)(1.1)(8.8) 


1.15(09)(14) 
0.83(09)(12) 
0.21(09)(13) 
0.47(08)(11) 


6.47(51) (78) (04) x 10 3 
4.64(53)(69)(03) x 10 3 
1.19(48)(71)(08) x 10 3 
2.62(44)(60)(02) x 10 3 



volume dependence. As the nucleon is comprised of light valence quarks only, it is expected 
to couple most strongly to pions, which should dominate its finite-volume modifications in 
the large- volume limit. It is expected that baryons with more strange quarks exhibit less 
volume sensitivity, and that the finite-volume mass shifts to the baryons, 5(FV)b, should 
naively satisfy the hierarchy 



S(FV) N > S(FV) E , 5(FV) A > 6{FV) S 



(15) 



II 



are shown in fig. 



The fit coefficients, c^, given in table 
approximately observed within the uncertainties of the ca 
of the E is somewhat smaller than naive expectations would suggest. 



where the expected hierarchy is 
culation. The volume dependence 
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0.211 
0.210 
0.209 
z 0.208 
-°~ 0.207 
0.206 
0.205 
0.204 

0.233 
0.232 
0.231 
0.230 

w 

S 0.229 
xf 

0.228 
0.227 
0.226 
0.225, 



♦ M N [m^ ~ 390 MeV] 

la stat 

la stat+syst 



0.001 0.002 0.003 0.004 0.005 0.006 

Expt-m^L] / m n h 



♦ M E [m^ ~ 390 MeV] 

la stat 

la stat+syst 




0.227 
0.226 
0.225 
0.224 
0.223 
0.222 
0.221 
0.220 
0.219, 



♦ M A [m n ~ 390 MeV] 

la stat 

la stat+syst 



0.001 0.002 0.003 0.004 

Expl-m^L] / m„L 



0.005 0.006 



0.001 0.002 0.003 0.004 0.005 0.006 

Expl-m^L] / rrijjL 



0.246 
0.245 
0.244 
0.243 - 

[■] 

2 0.242 - 
xf 

0.241 - 
0.240 
0.239 
0.238, 



♦ M s [m^ 390 MeV] 

- la stat 
la stat+syst 




0.001 0.002 0.003 0.004 

Expt-m^L] / rrijjL 



0.005 0.006 



FIG. 2: The mass of the nucleon (upper left panel), the A (upper right panel), the £ (lower left 
panel) and the H (lower right panel) as a function of e~ m ^ L / {m^L). The points and associated 
uncertainties (blue) are the results of the Lattice QCD calculations. The dark (light) shaded region 
corresponds to the la statistical uncertainty (statistical and systematic uncertainties combined in 



quadrature) associated with a fit of the form given in eq. (14). 



B. The Octet Baryons with NLO HB%PT 

In this section we explore both SU(2) L ® SU(2) R and SU(3) L ® SU(3) H HB^PT predictions 
and fits to the results of the Lattice QCD calculations. The analyses are performed at NLO 
in the chiral expansion; unfortunately they do not provide significant constraints on the 
counterterms that appear beyond NLO in HB^PT. Our strategy in these analyses is to use 
the octet-octet axial couplings and the octet-decuplet mass splittings from experimental data 
and Lattice QCD results, and fit the octet-decuplet axial couplings and the baryon masses 
in the infinite-volume limit to the results of the Lattice QCD calculations, given in table 
using two-flavor HB^PT. Inserting these fit values into the full three-flavor finite-volume 
corrections gives a measure of the relevance of kaon and rj loops. The goal is to determine 
the extent to which two- and three-flavor HB;\PT describe the volume dependence of the 
results of the Lattice QCD calculations. And, of course, it is of interest to determine whether 
any significant constraints can be placed on the -with few exceptions, rather poorly known- 
axial coupling constants of the baryons by studying finite- volume effects. 
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FIG. 3: The fit values of 4 , the coefficient of e - m ' rL /(m vr L) in eq. dl#, given in table 
each of the octet baryons. 



II 



for 



1. Parameter Set 

Lattice QCD calculations show that the nucleon axial coupling, g^, is essentially indepen- 
dent of the light-quark masses [28H3T] . and so in the following we use the experimentally- 
determined value qa = 1-26 as well as the central values of #ae, <?ee, and gss determined 
from Lattice QCD calculations [32] interpolated to the appropriate pion mass, or from the 
tree-level SU(3) relations between axial couplings at the physical pion mass (which are con- 
sistent with each other) 

p AS = 1.58(20) , <fe E = 0.900(30) , g 3S = 0.262(13) . (16) 

While the pion decay constant, f„, is experimentally determined to be f n ~ 132 MeV at 
the physical light-quark masses, Lattice QCD calculations have determined how it depends 
upon the light-quark masses, and at m w ~ 390 MeV its value is f w ~ 150 MeV • We 

take the baryon mass splittings determined from octet and decuplet correlation functions 
calculated on the 32 3 x 256 ensemble 4 : 

A AN = 298 MeV , A EN = 128 MeV , A AN = 90 MeV , A S , N = 427 MeV ; 
A SA = 38 MeV , A S , A = 336 MeV , A SA = 108 MeV , A S » A = 406 MeV ; 
A HS = 69 MeV , A s « s = 298 MeV , A as = 229 MeV , A H *s = 368 MeV ; 
A S * H = 298 MeV , A S * H = 228 MeV , A n3 = 367 MeV . (17) 

For the SU(3) L <g> SU(3) R HB^PT analysis, the SU(3)-symmetric axial couplings are 
fixed to the central values of the best-fit experimental values: D = 0.79, F = 0.47 and 



4 We do not quote uncertainties on these determinations as they do not significantly affect the NLO HB^PT 
fits. 
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C = 1.47 [37]. Further, the decay constants and masses are set to fK = f v = 160 MeV and 
the lattice-determined values = 544 MeV, and m v = 587 MeV (the latter determined 
from m n and nix via the Gell-Mann-Okubo (GMO) relation), respectively. 



2. The Nucleon Mass 
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FIG. 4: The Nucleon EMP's determined on the four lattice ensembles used in this work. They 
each result from linear combinations of different correlation functions that optimize the plateau of 
the ground state. Note that the y-axis scale is the same in all four panels. 



The nucleon EMP's and fits to the mass plateaus obtained from the results of the Lattice 
QCD calculations on the four lattice ensembles are shown in fig. |4j With the input pa- 
rameters given in subsection |IVB 1[ a two-parameter fit to the nucleon mass data, given in 



table [II can be performed in two-flavor HBxPT to determine the infinite-volume value of 
the nucleon mass, , and the value of the AiV7r axial coupling, g^Nn- While there have 
been many previous determinations of gANn, we treat it as a fit parameter and compare its 
value with the previous extractions. Fitting the expression in eq. ^ to the results given in 
table H gives the fit-regions shown in fig. [5j and fit values of 



M 



(oo) 
N 



0.20455(19)(17) t.l.u , I^a^I = 2.80(18)(21) . 



(18) 
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The ratio of couplings 

Ifl'ATVTrl 



9 a 



2.22(14)(17) , (19) 



can be compared with the value of \gANn\/9A = 1.56(06) [38] 5 extracted from an analysis of 
experimentally-measured nN scattering phase shifts. The difference between these two val- 
ues is significant, but as the two extractions have been performed at two different unphysical 
pion masses, little can be concluded. On the other hand, direct Lattice QCD calculations of 
\9antt\/9a have been performed [311110] over a range of pion masses 6 . One such calculation 
performed with a pion mass in the vicinity of ~ 390 MeV gives \gANn\/gA = 1.47(19). 
The finite-volume corrections to the nucleon mass resulting from this value of the coupling is 
shown as the dashed (red) curve in fig. [5j and clearly the contribution from the Nn and the 
An intermediate states constitute a large fraction of the finite-volume shift of the nucleon 
mass. Given the size of the pion mass in the present calculations, m n /M^ ~ 0.35, we antic- 
ipate that higher orders in HB^PT will change the finite-volume corrections at the ~ 30% 
level, consistent with the difference between the results of the lattice QCD calculations and 
NLO in HB^PT 7 . A next-to-next-to-leading order (NNLO) calculation of the finite-volume 
contributions to the nucleon mass in HB^PT, accompanied by more precise Lattice QCD 
calculations over a range of lattice volumes and quark masses, is required in order to improve 
upon this determination of gANn 8 - 



With these fit parameters and the parameter set previously defined in subsection |IVB 1 
the effect of kaon and r\ loops can be estimated by including the finite-volume corrections 
given in eq. the results of which are shown in fig. [5] The contributions from the strange- 
baryon and strange-meson intermediate states are estimated to be small, and somewhat 
improve the agreement between theory and the Lattice QCD calculation. Including them in 
the fit of g^Nn to the results of the Lattice QCD calculation gives \gANn\/9A = 2.10(15)(20), 



which is to be compared with the result in eq. (19). 

As NLO HBxPT reproduces most of the volume dependence of the nucleon mass at 
~ 390 MeV, and is expected to become more accurate at lighter pion masses, it is 
useful to use the NLO expression to estimate the size of the finite-volume contributions 
to the nucleon mass at the pion masses other than the current one. In fig. [6] we show 
the finite-volume contributions to the nucleon mass that are predicted at NLO in HB^PT 
for ~ 390, 230 and 140 MeV. This makes clear that finite-volume effects are expected 



5 The value of \gANn\/9A m R-ef. |38j has been divided by v2 in order to match the definition of gANv 
employed in defining eq. (p| |17j . 

6 The two methods employea in Ref. [40] suggest that there may be a relatively large systematic uncertainty 
in their value of \gANir\/gA beyond that quoted. 

7 If instead of using f v ~ 150 MeV to evaluate the NLO HBxPT result, the value at the physical pion mass, 
/,r ~ 132 MeV is used, then a value of \gANn\l <M — 1-76(18) is obtained, consistent within uncertainties 
with the extraction from the matrix element of the axial current |39| . The ambiguity in the value of 
the decay constant that is used in the NLO contribution will be parametrically reduced by a NNLO 
calculation. 

8 An alternative, efficient way of capturing the bulk of the volume dependence is to insert the forward 7rN 
scattering amplitude, calculated in HBxPT with explicit A degrees of freedom into eq. [41J . This 
procedure has been shown to work well for the n mass volume dependence |42j . 
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FIG. 5: The mass of the nucleoli as a function of e m " ri /(m 7r L). The dark (light) grey shaded re- 
gion corresponds to the la statistical uncertainty (statistical and systematic uncertainties combined 
in quadrature) resulting from fitting and gANw Using this value of M^°\ the dot-dashed 

curve (green) corresponds to the contribution from octet baryons and pions, the dotted curve (or- 
ange) corresponds to the contribution from octet baryons and kaons or an rj, the dashed curve 
(red) corresponds to the contribution from octet and decuplet baryons and pions, and the solid 
curve (blue) corresponds to the contribution from octet and decuplet baryons and pions, kaons or 
an rj. 



to be significantly smaller at the lighter pion masses for fixed m^L. The values of m T L, 
estimated at NLO in HB;\PT, for which the finite-volume contributions to the nucleon 
mass are 5M ( ^ V) = 1 MeV at m n = 390,230 and 140 MeV are m n L ~ 6.2,4.7 and 3.9, 

respectively. For 8M^ = 100 keV, the corresponding values are m^L ~ 8.0, 6.4, and 
5.8, respectively. Given that the deuteron binding energy and nuclear excitation energies 
are in the MeV regime, the estimates indicate that m„-L> 2n is required at the physical 
pion mass in order to eliminate contamination from this class of exponentially-suppressed 
finite-volume effects that contaminate the extraction of phase shifts and binding energies 
from Lattice QCD calculations. 



3. The A Mass 

The A EMP's and fits to the mass plateaus obtained from the results of the Lattice QCD 
calculations for the four lattice ensembles are shown in fig. [7j The fit values of the A masses 



14 




FIG. 6: Estimates of the finite- volume contributions to the nucleon mass at NLO in HB^PT 
for rn-Tr = 390 MeV (upper, blue, solid, which corresponds to the red dashed curve in fig. [5]), 
m n = 230 MeV (middle, green, dotted), and m n = 140 MeV (lower, red, dashed). 



in the four lattice volumes are given in tablejl], and are shown as the points with uncertainties 
in fig. [8] The shaded regions in fig. [8] show the results of the SU(2) L ® SU(2) R HB^PT fit to 
the volume dependence of the A mass using eq. ^ and the value of g^A given in eq. (16). 
The fit gives 

M| oo) = 0.22064(15)(19) t.l.u , \g^ Aw \ = 2.21(16)(23) . (20) 

If instead of fitting gs*Awi flavor SU(3) symmetry is used to relate it to gANn, 9t,*a-k = 
9anit/V2 — 1-3, then the contribution from S intermediate states and from S and E* 
intermediate states are shown as the dot-dashed (green) and dashed (red) curves in fig. [8j 
respectively. Comparing these expectations with the results of the Lattice QCD calculations, 
manifested in the fit value of g^An being ~ 40% larger than phenomenological expectations, 
indicates that higher orders in two-flavor %PT are important, or that the strange quark 
plays a role in the finite- volume contributions through kaons or an rj. 

As with the nucleon, we can now estimate the effects of kaon and r\ loops by including 



the finite- volume corrections given in eq. (10). This gives rise to the curves shown in fig. 8 



(Note that the same values for the input parameters are used to generate the curves in 
fig. [8] as are used in the case of the nucleon) . The dot-dashed (green) and dashed (red) 
curves are the contributions from pions, while the dotted (orange) and solid (blue) curves 
correspond to the sum of contributions from pions, kaons and an rj. The pions make the 
largest contribution to the volume dependence of the mass of the A, but, unlike in the case 
of the nucleon, the kaons and rj contribute significantly. It is interesting to note that the 
NLO contribution in three-flavor HBxPT (blue curve) agrees well with the results of the 
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FIG. 7: The A EMP's determined on the four lattice ensembles used in this work. They each result 
from linear combinations of different correlation functions that optimize the plateau of the ground 
state. Note that the y-axis scale is the same in all four panels. 

Lattice QCD calculation. However, NLO HB%PT is expected to be modified at the ~ 30% 
level by higher orders in the expansion. 



4- The £ Mass 

The £ EMP's and fits to the mass plateaus obtained from the results of the Lattice QCD 
calculations on the four lattice ensembles are shown in fig. |9j The fit values of the S masses 
in the four lattice volumes are given in table [TJ, and are shown as the points with uncertainties 
in fig. 10 The shaded regions in fig. 10 show the results of the SU(2) L eg) SU(2) H HB^PT 
fit to the volume dependence of the £ mass using eq. ^ with the axial couplings given in 
eq. (16). The fit gives 



M 



(oo) 



0.22747(17) (19) t.l.u 



\g^\ < 1.38(1.90] 



(21) 



where we have quoted a 68% confidence interval for ^e'Stt including statistical errors and 
statistical and systematic errors added in quadrature (bracketed). If instead of fitting g£*£7r, 
the SU(3) relation is used, g^^n = gANw/V3 = 1.07, then the contribution from E, A 
intermediate states and from S, A and £* intermediate states are shown as the dot-dashed 
(green) and dashed (red) curves in fig. flOj respectively. The NLO SU(3) L £g> SU(3) H HB^PT 
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FIG. 8: The mass of the A as a function of e m " ri /(m 7r L). The dark (light) grey shaded region 
corresponds to the la statistical uncertainty (statistical and systematic uncertainties combined in 
quadrature) resulting from fitting Ml and <7e*Att- Using this value of m|°° , the dot-dashed curve 
(green) corresponds to the contribution from octet baryons and pions, the dotted curve (orange) 
corresponds to the contribution from octet baryons and kaons or an rj, the dashed curve (red) 
corresponds to the contribution from octet and decuplet baryons and pions, and the solid curve 
(blue) corresponds to the contribution from octet and decuplet baryons and pions, kaons or an 77. 



prediction (blue curve) is somewhat higher than the result of the Lattice QCD calculation 
in the smallest volume, but not significantly so. It should be added that the SU(3) value 
of the coupling constant, gs*s n = 1-07 is consistent with the confidence interval extracted 



from the SU(2) L ® SU(2) H fit, given in eq. Q. 

The volume dependence of the £ mass is found to be somewhat smaller than that of the 
A mass. This is consistent with the prediction of NLO SU(3) L (g> SU(3) R HB%PT, which 
is largely driven by the coupling to the decuplet intermediate states. The difference in the 



couplings to the decuplet, given in eq. (13), is sufficient to explain the difference in volume 
dependence. 



5. The S Mass 

The H EMP's and fits to the results of the Lattice QCD calculations on the four lattice 



ensembles are shown in fig. 11 The fit values of the S masses in the four lattice volumes are 
given in table [T| and are shown as the points with uncertainties in fig. 12 The shaded regions 
in fig. 12 show the results of the SU(2) L <g> SU(2) fl HB^PT fit to the volume dependence of 
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FIG. 9: The £ EMP's determined on the four lattice ensembles used in this work. They each 
result from linear combinations of different correlation functions that optimize the plateau of the 
ground state. Note that the y-axis scale is the same in all four panels. 



the S mass using eq. (|8j). The fit parameters are 
r(oo) = .23978(12)(18) t.l.u , 



Mi 



2.49(23)(35) 



(22) 



If instead of fitting gz*E.n, the estimates from other observables are used, gs*Sw = 9an-k/2 = 
0.93 [43J, the contribution from 5 intermediate states and from S and 5* intermediate states 



are shown as the dot-dashed (green) and dashed (red) curves in fig. 12 , respectively. Com- 
paring the expectations with the results of the Lattice QCD calculations, manifested in the 
fit value of <7s*Htt being more than twice expectations, suggest that the pionic contributions 
do not dominate the finite-volume corrections, even after considering contributions from 
higher orders in SU(2) L £g> SU(2) R HBxPT. As the H carries two strange quarks, one expects 
kaon and rj loops to make relatively larger finite-volume contributions to the S mass than 
to the nucleon, A and £ masses. 



The NLO SU(3) L ®SU(3) J? HB^PT prediction from eq. (12) gives rise to the curves shown 
in fig. 12 The dot-dashed (green) and dashed (red) curves are the contributions from pions, 
while the dotted (orange) and solid (blue) curves correspond to the sum of contributions 
from pions, kaons and an r\. HB%PT predicts that it is the kaons and r\ that dominate the 
finite-volume contributions to the H mass. The exponential suppression of the kaon and 
rj contributions is not sufficient to overcome the relatively large axial coupling constants 
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FIG. 10: The mass of the E as a function of e _m " rL /(m 7r L). The dark (light) grey shaded region 
corresponds to the la statistical uncertainty (statistical and systematic uncertainties combined in 
quadrature) resulting from fitting M^°°^ and 5s*Stt- Using this value of M^°\ the dot-dashed curve 
(green) corresponds to the contribution from octet baryons and pions, the dotted curve (orange) 
corresponds to the contribution from octet baryons and kaons or an rj, the dashed curve (red) 
corresponds to the contribution from octet and decuplet baryons and pions, and the solid curve 
(blue) corresponds to the contribution from octet and decuplet baryons and pions, kaons or an rj. 

at this pion mass. These results suggest that NLO SU(3) L Cg> 811(3)^ HB^PT provides a 
good description of the finite-volume modifications to the 5 mass. Again, one expects NLO 
HB^PT to provide an estimate of the finite volume effects that would be modified at the 
~ 30% level by higher orders in the expansion. 



6. Summary 



The infinite-volume values of the octet-baryon masses fit from the Lattice QCD data using 
HB^PT are consistent (almost identical) with the infinite-volume values extracted from 
the simple phenomenological fits presented above (see table III for a summary). However, 
there are several important lessons that one learns from the HB^PT analysis of the volume 
dependence of the octet-baryon masses. First, by comparing, for instance, the dashed (red) 
and dot-dashed (green) curves in fig. [5j fig. |8j fig. 10 and fig. 12, one sees the relevance 
of the octet-decuplet axial transitions in the description of the finite-volume effects. We 
conclude that HBxPT with the decuplet states integrated out does not provide a reliable 
description of the finite-volume effects at NLO at this pion mass. Second, by comparing, 
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FIG. 11: The 5 EMP's determined on the four lattice ensembles used in this work. They each 
result from linear combinations of different correlation functions that optimize the plateau of the 
ground state. Note that the y-axis scale is the same in all four panels. 



for instance, the solid (blue) curve and the dashed (red) curve in fig. [5j fig. [8j fig. 10 and 
fig. [12J one sees the relative importance of fluctuations to intermediate states involving kaons 
and/or 77, which are not captured in the two- flavor chiral expansion. While these effects are 
small in the case of the nucleon, they are significant for the hyperons at the heavy pion 
mass at which the Lattice calculations were performed. Therefore, while SU(2) L <g> SXJ(2) R 
HB^PT is adequate for the nucleon, SU(3) L ® SU(3) H HB%PT is necessary to account for 
the finite- volume effects experienced by the hyperons 9 . And thus, the fit hyperon axial 



couplings presented in table [Hl| should not be considered reliable. This is, of course, due to 
rriK — rn-K not being small enough, and therefore the SU(2) L £g> 811(2)^ fits for the hyperons 
will improve as the physical pion mass is approached. It is important to emphasize that we 
have not propagated the uncertainties associated with the input axial couplings and baryon 
mass splittings, as we have found that the resulting uncertainties in the fit quantities are 
smaller than the expected size of NNLO effects. The various curves in fig. [5j fig. [8j fig. 10 
and fig. [12] become bands when the input parameter uncertainties are included. 



9 We caution that reproducing the calculated volume dependence using HB^PT does not imply that the 
quark-mass dependence itself is under control. Indeed it is well known that SU^)^ ® SU(3) fl HB^PT 
converges poorly for certain quantities [TIH46| . 
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FIG. 12: The mass of the 5 as a function of e m " rZ '/(m 7r L). The dark (light) grey shaded region 
corresponds to the la statistical uncertainty (statistical and systematic uncertainties combined in 
quadrature) resulting from fitting Mi 00 ^ and 5h*Htt- Using this value of m1°°\ the dot-dashed curve 
(green) corresponds to the contribution from octet baryons and pions, the dotted curve (orange) 
corresponds to the contribution from octet baryons and kaons or an rj, the dashed curve (red) 
corresponds to the contribution from octet and decuplet baryons and pions, and the solid curve 
(blue) corresponds to the contribution from octet and decuplet baryons and pions, kaons or an rj. 



TABLE III: The results of NLO SU(2) L <g) SU(2) R HBxPT fits to the volume dependence of the 
baryon masses. M^°^ is the infinite- volume extrapolation of the baryon mass. The first uncertainty 
is statistical, the second is the fitting systematic, and the third (where appropriate) is due to scale 
setting. The uncertainties in the input axial couplings and baryon mass splittings have not been 
propagated, as discussed in the text. 



Hadron 


M^° ] (t.l.u.) 


M { ™ ] (MeV) 


Axial Coupling 


M N 


0.20455(19)(17) 


1151.3(1. 1)(1.0)(7.5) 


IsajvttI = 2.80(18)(21) 


M A 


0.22064(15) (19) 


1241.9(0.8)(1.1)(8.1) 


\gx*An\ = 2.21(16)(23) 


M E 


0.22747(17) (19) 


1280.3(1.0)(1.1)(8.3) 


ISE'&rl < 1.38[1.90] 


M s 


0.23978(12)(18) 


1349.6(0.7)(1.0)(8.8) 


Iph-HttI = 2.49(23)(35) 
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C. Combinations of Masses 

In addition to examining the volume dependence of the baryon masses, it is interesting to 
explore the volume dependence of certain combinations of the masses. In order to minimize 
both statistical and systematic uncertainties in determining various combinations of masses 
from the Lattice QCD calculation, in each case a correlation function is formed from prod- 
ucts or ratios of the individual baryon correlation functions under a jackknife or bootstrap 
procedure, from which the combination of masses is extracted. 



1. The Centroid of the Octet 

The centroid of the baryon octet is the sum of the masses weighted by the isospin degeneracy 
of each state, 

M 8 = ^M A + -M s + -M N + . (23) 



The results of the Lattice QCD calculations are shown in fig. 13, along with the results of 



a simple fit of the form described in subsection |IV A[ The simple fit gives rise to a centroid 
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mass of M 8 (oo) = 0.22354(17)(20) t.l.u = 1255.1(1.0)(1.1)(8.2) MeV. Also shown in fig. _ 
are the predictions of both SU(2) L <g> SU(2) fl and SU(3) L <g> SU(3) fl NLO HB^PT resulting 
from the same input parameters (not the fit parameters) as those used for the predictions 
of the individual baryon masses. In particular, the SU(3) relations between the decup let- 
octet axial coupling constants, and between the octet-octet kaon and rj axial couplings, 
are employed. Given the overall general agreement between the leading predictions and 
the individual baryon masses, it is unsurprising that the SU(3) L <S> SU(3) R prediction for 
the centroid mass of the octet agrees reasonably well with the results of the Lattice QCD 
calculation. 



2. The E-A Mass Splitting 

The mass difference between the E and the A vanishes in the limit of exact SU(3) symmetry. 
In Nature, the splitting is found to be M s — Ma ~ 74 MeV, and consequently, at the pion 
mass of the Lattice QCD calculations, the calculated splitting is expected to be small. The 



results of the Lattice QCD calculation are shown in fig. 14, along with the results of a 



simple fit, of the form given in eq. (14), shown as the shaded regions. The result of the 
simple fit gives (M s - M A ) (oo) = 0.006598(48)(63) t.l.u = 37.05(27)(35)(24) MeV, which is 
approximately half of its value at the physical quark masses. This is, in part, due to the 
strange quark mass used in the calculation being somewhat lighter than its value in Nature. 
The finite-volume contributions significantly suppress the mass splitting in smaller volumes. 

The NLO expressions for the mass splitting do not describe the observed volume depen- 
dence well. While the full SU(3) L ® SU(3) R NLO amplitude agrees in its sign with the 
volume dependence, the magnitude is significantly smaller than the results of the Lattice 
QCD calculations. It is clear that SU(3) breaking contributions that enter at higher orders 
in the chiral expansion play an important role in the volume dependence of the E-A mass 
splitting. 



22 



0.228 



0.227 



0.226 



0.225 



0.224 



0.223 



0.222 



Im^- 390 MeV] 



♦ M, 

8 baryons + K/r| 

— 8 and 10 baryons + Jt + K/n 

— 8 and 10 baryons + % 
■ — 8 baryons + n 

la stat 

la stat+syst 




_L 



0.001 0.002 0.003 0.004 

Expf-nijjL] / rn^L 



0.005 



0.006 



FIG. 13: The centroid of the baryon octet as a function of e~ m7rL /(m^L). The dark (light) 
grey shaded region corresponds to the la statistical uncertainty (statistical and systematic uncer- 
tainties combined in quadrature) resulting from fitting the m n L = oo value and the coefficient of 
e -m n L J (jn^^y Using the value at m^L = oo, the dot-dashed curve (green) corresponds to the con- 
tribution from octet baryons and pions, the dotted curve (orange) corresponds to the contribution 
from octet baryons and kaons or an r], the dashed curve (red) corresponds to the contribution from 
octet and decuplet baryons and pions, and the solid curve (blue) corresponds to the contribution 
from octet and decuplet baryons and pions, kaons or an rj. 



3. The Gell-Mann-Okubo Mass Relation 
The GMO relation, 

GMO = M A + ' r \/v - 2 -M N - 2 -M R , (24) 

vanishes in the limit of exact SU(3) flavor symmetry, and also vanishes in the limit where 
the SU(3) breaking transforms as an 8 under SU(3) flavor symmetry. Consequently, it 
is a valuable probe of the structure of flavor symmetry breaking, being non-zero only for 
breaking that transform in the 27 irreducible representations of SU(3) 10 . The results of the 



Lattice QCD calculations are shown in fig. 15, along with the results of the simple fit, of 



the form given in eq. (14), shown as the shaded regions. The result of the simple fit gives 



10 Only the symmetric irreps in 8 <g> 8 = 27 © 10 © 10 © 8 © 8 © 1 are allowed, i.e. the 27, 8, and 1. By 
construction the 8 and 1 do not contribute. 
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FIG. 14: The E-A mass difference as a function of e~ m ^ L j (rn^V) . The dark (light) grey shaded 
region corresponds to the 1<7 statistical uncertainty (statistical and systematic uncertainties com- 
bined in quadrature) resulting from fitting (Mj — Ma)^ 00 -* and the coefficient of e~ mnL /(m-^L). 
Using this value of (Ms — M\)^°°\ the dot-dashed curve (green) corresponds to the contribution 
from octet baryons and pions, the dotted curve (orange) corresponds to the contribution from octet 
baryons and kaons or an 77, the dashed curve (red) corresponds to the contribution from octet and 
decuplet baryons and pions, and the solid curve (blue) corresponds to the contribution from octet 
and decuplet baryons and pions, kaons or an r/. 



GMO (oo) = 3.49(25)(46) x 1(T 4 t.l.u = 1.96(14)(26)(01) MeV. Given the smallness of the 
GMO combination, it is no surprise that it has a large volume dependence, even changing 
sign between the 20 3 x 128 and the 24 3 x 128 lattice volumes. 

As the GMO relation is sensitive only to the 27 SU(3) breaking, both the finite-volume 
and infinite- volume contributions are calculable in SU(3) L <S> SU(3) R HBxPT at one loop 
(which gives a finite result). It is straightforward to show that the infinite- volume value of 
the GMO relation is H9l 



GMO (NLO) 



247T/2 

where the function F c = F(m c , A, //) is 



Ql>2 _ (4m 3 , - 3m 3 - m 3 ) - ^ {AF K - 3F V - F w 



■(25) 



F(m,A,/i) = (m 2 -A 2 ) fy^^log E + ie ) - A log ( ^ 



/< 2 



\^\og{^) . (26) 
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FIG. 15: The GMO relation as a function of e _m ' ri /(m 7r L). The dark (light) grey shaded region 
corresponds to the la statistical uncertainty (statistical and systematic uncertainties combined in 
quadrature) resulting from fitting the m^L = oo value and the coefficient of e~ m7rL /(m n L). The 
dot-dashed curve (green) corresponds to the contribution from octet baryons and pions, the dotted 
curve (orange) corresponds to the contribution from octet baryons and kaons or an 77, the dashed 
curve (red) corresponds to the contribution from octet and decuplet baryons and pions, and the 
solid curve (blue) corresponds to the contribution from octet and decuplet baryons and pions, 
kaons or an yj. 



We have set f-w = fx = 150 MeV in the GMO relation to eliminate formally higher-order 
contributions that depend upon the renormalization scale, /1, from this expression. Inserting 



the values of the constants and parameters that we have used previously into eq. (26) gives a 
value of GMO (iVLO) ~ 6.8 x 10" 4 t.l.u., which is a factor of two greater than the extrapolation 
of the Lattice QCD results. This is not that surprising given the expected size of the 
higher-order contributions in the HB%PT expansion, as well as the fact that this quantity 
is anomalously small (suppressed beyond naive expectations based upon SU(3) symmetry 
alone due to a further suppression by l/N% in the large- iV c limit of QCD n ). 

The NLO predictions of the volume dependence of the GMO relation in SU(2) L <g> SU^)^ 
and SU(3) L <g) SU(3) fl HB^PT are shown in fig. 15 The contribution from the kaons and 



the r\ are of opposite sign to that from the pion (as expected from their cancellation in the 



SU(3) limit), which gives rise to a small volume dependence, even on the scale of fig. 15 



m Ut d denotes the SU(3) breaking parameter, then the GMO relation scales as ~ e 3 / 2 /AT 



If e ~ n 

relative to the baryon masses; see, for example, eq. (25). For a review, see Ref. |T7] . 
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The predictions for the volume dependence of the GMO relation are in clear contradiction 
with the results of the Lattice QCD calculations. Clearly NLO HB^PT does not describe 
the higher-dimensional SU(3) breaking that provides the finite- volume dependence of the 
GMO relation, and we have found that this relation is particularly sensitive to the volume 
of the lattice. 

The GMO relation was previously explored by some of the present authors [IH HEJ • At 
approximately the same pion mass, and in the volume with L ~ 2.5 fm, the GMO relation 
was found to be positive and consistent with the loop-level expression given in eq. (26). 
This is in slight disagreement with the current determination; however, neither calculation 
extracts values in the continuum limit and the difference in this very small quantity may arise 
from the different discretizations and/or the somewhat different values of the strange-quark 
mass. 



4- The i?4 Mass Relation 
The i?4 relation, defined to be 

Ra = g ( M N + M A + M s - 3M S ) , (27) 

vanishes in the limit of exact SU(3) flavor symmetry, and is formally dominated by a single 
insertion of the light-quark mass matrix in the infinite volume limit. This relation is of phe- 
nomeno logical interest because not only does it vanish in the limit of exact SU(3) flavor sym- 
metry, but it also vanishes in the large- N c limit of QCD |47J SHJ ED] , scaling as ~ e/N c relative 



to the baryon masses (e is defined in the footnote in subsection IVC 3). Recently, this rela- 
tion, along with other relations among masses were examined with Lattice QCD calculations 
in work by Jenkins et al. [50] using domain-wall light-quark and strange-quark propagators 



generated on a number of ensembles of improved Kogut-Susskind dynamical quarks [5T| 152] . 



The results of the present Lattice QCD calculations are shown in fig. [16j along with the 
results of the simple fit, of the form given in eq. flli] ), shown as the shaded regions. The 
result of the simple fit gives = -0.002757(27) (40) t.l.u = -15.48(15) (22) (09) MeV, 

which is approximately half of its value at the physical quark masses, 

^expt = _ 34 5 MeV _ 

The finite-volume contributions significantly suppress the mass splitting in smaller volumes. 

The NLO expressions for R4 do not describe the observed volume dependence well. While 
the full SU(3) L ® 811(3)^ NLO amplitude agrees in its sign with the volume dependence, the 
magnitude is significantly smaller than the results of the Lattice QCD calculations. SU(3) 
breaking contributions that enter beyond NLO in HBxPT play an important role in the R4 
mass relation. 



V. THE VOLUME DEPENDENCE OF THE MESON MASSES 

It is also interesting to explore the volume dependence of the meson masses, specifically 
that of the pion and the kaon. There has been significantly more theoretical and numerical 
exploration of the meson masses and how they depend upon the volume of the lattices used 
in Lattice QCD calculations. An overview of the theoretical status of such finite-volume 
contributions can be found in Ref. [53J. 
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FIG. 16: The mass relation as a function of e~ m7lL / (m n L). The dark (light) grey shaded region 
corresponds to the la statistical uncertainty (statistical and systematic uncertainties combined in 
quadrature) resulting from fitting the m^L = oo value and the coefficient of e~ m7rL /(m n L). The 
dot-dashed curve (green) corresponds to the contribution from octet baryons and pions, the dotted 
curve (orange) corresponds to the contribution from octet baryons and kaons or an 77, the dashed 
curve (red) corresponds to the contribution from octet and decuplet baryons and pions, and the 
solid curve (blue) corresponds to the contribution from octet and decuplet baryons and pions, 
kaons or an 77. 



The results of the present Lattice QCD calculations of the meson masses in the four 



different lattice volumes are given in table |IV| and the EMP's are shown in fig. [17] and 
fig- 



TABLE IV: Meson masses from the Lattice QCD calculations in the four lattice volumes. 



L 3 x T 


16 3 x 128 


20 3 x 128 


24 3 x 128 


32 3 x 256 


m n (t.l.u.) 


0.06943(36) (0) 


0.06936(12)(0) 


0.06903(19)(0) 


0.069060(66)(81) 


rriK (t.l.u.) 


0.09722(26)(0) 


0.09702(10)(03) 


0.09684(15)(01) 


0.096984(78) (60) 
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FIG. 17: The pion EMP's determined on the four lattice ensembles used in this work. Note that 
the y-axis scale is the same in all four panels. 



A. The Pion Mass 



The finite-volume contribution to the mass of the pion in SU(2) L (g> SU(2) fl %PT is given 
bydU 



ra^{L) — m n (oo) 



3v 3 



(28) 



where K\(x) is the modified Bessel function. The meson masses have different overall volume 
scaling to the baryons, due to the absence of a three-meson vertex. As Kx(z) —> e~ z / y/z, the 



results of the Lattice QCD calculations are shown in fig. 19 as a function of e~ m ' ri /(m vr L) 3 / 2 
rather than e~ m " rL /(m vr L) as was used for the baryons. Consequently, the naive fit that we 
perform to the meson masses is of the form 

mZ\m n L) = mffi + ^ ^ 3/2 • (29) 

With the current precision of the Lattice QCD calculation, we cannot distinguish be- 
tween the fit forms of e _m " rL /(m 7r L) and e~ m ' ri /(m 7r L) 3 / 2 with statistical significance. 
The fit parameters are mi° o) = 0.069073(63)(62) t.l.u. = 387.8(0.4)(0.4)(2.5) MeV and 
c[ y ) = 0.23(12)(07) t.l.u. = (1.30(65)(39)(01)) x 10 3 MeV. 
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FIG. 18: The kaon EMP's determined on the four lattice ensembles used in this work. Note that 
the y-axis scale is the same in all four panels. 



At NLO in SU(3) L <8> 811(3)^ xPT, the finite-volume corrections to the pion mass are 
given by [53J 



5m m 



3ml 



,2 t 



m^m 



'/ 



3v 3 

Ki(m„L) + yplK^sPlm^V) + -^K x [^lm r] L) + ... 



.(30) 



47r 2 /^ rn v L 

Using the value of mi°°^ found in the fit to the form in eq. (|29l), the predicted volume 



dependence is shown in fig. 19 as the solid ( SU(3) i ®SU(3) J? ) and dashed curves ( SU(2) L £g> 
SU(2) R ). The volume dependence that is found in the Lattice QCD calculations agrees with 
the expectations of NLO xPT, and is significantly smaller than that of the baryon masses. 



B. The Kaon Mass 



The formalism describing the volume dependence of the kaon mass is analogous to that of 
the pion. The NLO calculation in SU(3) L <g> 811(3)^ xPT gives 

,2 i 



5m 



mj<m 



K 



'i 



2tt%2 mri L 



Ki(m v L) + y/2Kx(V2m n L) + —^K^VSm^L) + .. 

3y3 



.(31) 
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FIG. 19: The mass of the pion as a function of e " 1 ' ri /(m 7r L) 3 / 2 . The points and associated 



uncertainties (blue) are the results of the Lattice QCD calculations, as given in table IV The 



dark (light) shaded region corresponds to the la statistical uncertainty (statistical and systematic 



uncertainties combined in quadrature) associated with a fit of the form given in eq. ( 29 ) . The solid 
(red) curve corresponds to the prediction of NLO SU(3) L <8> SU(3) B xPT, while the dashed (green) 
curve corresponds to the prediction of NLO S\J(2) L ^ 1 1 ; ' '• >- r> ' r wv <>•'• " x 1 

in the fit. 



SU(2) R xPT using the value of m\ found 
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The results of the Lattice QCD calculatio n, g iven in table IV, are shown in fig. 
along with a fit to the form given in eq. (29). The resulting fit values are = 
0.096953(68)(39) t.l.u. = 544.4(0.4) (0.2) (3.5) MeV and c ( p = 0.087(93)(44) t.l.u. = 
(4.9(5.2) (2.5) (0.0)) x 10 2 MeV. The volume dependence of the kaon is found to be very 
small, much smaller than that of the baryons, and is consistent with the absence of pion 
loop contributions. 



VI. CONCLUSIONS 

We have performed precise Lattice QCD calculations of the low-lying hadron masses at a 
pion mass of m n ~ 390 MeV in four ensembles of anisotropic clover gauge-field configurations 
with a spatial lattice spacing of b s ~ 0.123 fm, an anisotropy of £ = 3.5 and cubic spatial 
lattice volumes with extent L ~ 2.0, 2.5, 3.0 and 3.9 fm. These calculations have allowed for 
a detailed exploration of the volume dependence of the octet baryon masses and of the pion 
and the kaon masses. 

Our main conclusions are as follows: 
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FIG. 20: The mass of the kaon as a function of e m ' ri /(?7i 7r L) 3 / 2 . The points and associated 



uncertainties (blue) are the results of the Lattice QCD calculations, as given in table IV The 



dark (light) shaded region corresponds to the la statistical uncertainty (statistical and systematic 

Df the f< 
SU(3) fi X PT. 



uncertainties combined in quadrature) associated with a fit of the form given in eq. ( 29 ) . The solid 
(red) curve corresponds to the prediction of NLO SU(3) L 



• In order to calculate individual baryon masses with percent-level precision (±10 MeV), 
it is sufficient to work in volumes with m^L^ 4.3 for ~ 390 MeV, and NLO HB^PT 
indicates that somewhat smaller values of m n L may be sufficient at lighter pion masses. 

• The expectations of NLO SU(3) L (g) SU(3) K HB^PT (with the meson decay constants 
evaluated at the appropriate meson mass) are found to qualitatively reproduce the 
lattice results for the volume dependence of all of the baryon and meson masses. The 
NLO SU(3) L (g) SU(3) R predictions are sufficiently accurate to allow for meaningful 
extrapolations to lighter pion masses to be made, where the volume dependences 
are expected to be significantly smaller. In the context of determinations of two- 
body interactions, this feature allows for gauge-field configurations with somewhat 
(logarithmically) smaller values of m n L to be used for the calculation of the interactions 
between baryons while keeping the exponential corrections to the Liischer eigenvalue 
relation from single hadron masses negligibly small. A Lattice QCD calculation at the 
physical pion mass and in a volume with m n L = 5.8 is predicted to generate finite- 
volume corrections to the nucleon mass of SMn ~ 100 keV, smaller than the typical 
nuclear excitation energies found in light nuclei. 

• The contributions from kaon and r\ loops to the volume dependence of the nucleon mass 
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are small in SU(3) i <E>SU(3) i j HB^PT, and therefore a relatively stable determination of 
\9antt\/9a has been found, which can be systematically improved by working at higher 
orders in HB%PT. Given the size of the finite-volume contributions to the masses of 
the hyperons from kaon and r] loops, we have not made significant determinations 
of the hyperon axial coupling constants in the two-flavor chiral expansion. However, 
Lattice QCD calculations at lighter pion masses should enable a determination of these 
couplings. 

• The GMO relation is found to exhibit substantial volume dependence, with the re- 
lation changing sign at m T L ~ 5.2 with ~ 390 MeV. Perhaps this should not 
be surprising given the delicate cancellations that occur between the baryon masses 
to leave a quantity that depends only upon SU(3)-breaking in the 27-dimensional ir- 
reducible representation and vanishes in the large- N c limit as 1/N% compared to the 
individual baryon masses. 
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